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Abstract. In this paper, we use the theory of symmetric Dirichlet forms to derive 
Feynman-Kac formulae for the forward problem of electrical impedance tomogra¬ 
phy with possibly anisotropic, merely measurable conductivities corresponding to 
different electrode models on bounded Lipschitz domains. Subsequently, we em¬ 
ploy these Feynman-Kac formulae to rigorously justify stochastic homogenization 
in the case of a stochastic boundary value problem arising from an inverse anomaly 
detection problem. Motivated by this theoretical result, we prove an estimate for 
the speed of convergence of the projected mean-square displacement of the under¬ 
lying process which may serve as the theoretical foundation for the development 
of new scalable stochastic numerical homogenization schemes. 


1. Introduction 

Electrical impedance tomography (EIT) aims to reconstruct the unknown conduc¬ 
tivity K in the conductivity equation 

(1) V • (kVu) = 0 in U 

from current and voltage measurements on the boundary of the domain D. This inverse 
conductivity problem is known to be severely ill-posed, that is, its solution is extremely 
sensitive with respect to measurement and modeling errors. As a result, EIT suffers 
from inherent low resolution and due to this limitation, many practical applications 
focus on the detection of conductivity anomalies in a known background conductiv¬ 
ity rather than conductivity imaging. In the mathematical modeling of such inverse 
anomaly detection problems, randomness typically reflects a lack of precise informa¬ 
tion about the meso- and microstructure of the heterogeneous background conductivity, 
which may fluctuate on many scales. Recently, the second author has proposed a novel 
method for the detection of conductivity anomalies in a random background conduc¬ 
tivity which is based on homogenization of the underlying stochastic boundary value 
problem, cf. [55]. 

Although the homogenization theory for elliptic divergence form operators is well- 
developed, cf., e.g., [6, 47, 48, 64], the numerical approximation of the effective con¬ 
ductivity in the random setting still poses major challenges. The commonly used 
deterministic methods based on a discretization of the so-called auxiliary problem have 
two main drawbacks. First, the auxiliary problem is formulated on the whole space 
and second, it has to be solved for almost every realization of the random medium. 
That is, truncations of the auxiliary problem have to be considered and choosing an 
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appropriate spatial truncation with appropriate boundary conditions is a delicate is¬ 
sue, cf. [11]. Moreover, in practically relevant cases, such as high contrast digitized 
random media, it is extremely difficult to solve the corresponding variational prob¬ 
lems by usual deterministic methods, such as the finite element or the finite difference 
method, due to the behavior of the solutions near the corner points. Therefore, numer¬ 
ical approximation of the effective conductivity can be prohibitively expensive in terms 
of computation time. As a matter of fact, practitioners often choose to avoid these 
computations at all and rather content themselves with theoretical bounds, cf., e.g., 
[60]. However, it has been reported in the physical literature that the shortcomings of 
the standard deterministic methods can be circumvented by using continuum micro¬ 
scale Monte Carlo simulation of certain diffusion processes evolving in random media 
instead, cf., e.g., [36, 56, 61, 38, 39]. In this work we give a rigorous mathematical jus¬ 
tification for homogenizing the EIT forward problem using such methods by studying 
the interconnection between reflecting diffusion processes and certain boundary value 
problems for the conductivity equation. More precisely, we derive Feynman-Kac for¬ 
mulae for solutions of the deterministic conductivity equation (1) posed on a bounded 
domain D C Mf", d > 2, with Lipschitz boundary dD and possibly anisotropic uniformly 
elliptic and uniformly bounded conductivity subject to different boundary conditions 
modeling electrode measurements. Subsequently, we employ the Feynman-Kac formu¬ 
lae to prove a homogenization result for the corresponding stochastic boundary value 
problem which justifies the use of stochastic numerical homogenization schemes based 
on simulation of the underlying diffusion processes in order to approximate the effec¬ 
tive conductivity. Finally, we prove an estimate for the speed of convergence of the 
projected mean-square displacement of the underlying diffusion processes. The main 
advantage of the presented approach to numerical homogenization, beside its inherent 
parallelism, is that its convergence rate is dimension-independent and its computational 
cost grows only linearly with the dimension. 

It is well known, that reflecting diffusion processes generated by non-diver-gence 
form operators with smooth coefficients on bounded, smooth domains, are Feller pro- 
eesses satisfying Skorohod type stochastic differential equations. The construction in 
the case of divergence form operators with merely measurable coefficients requires the 
theory of symmetric Diriehlet forms which has its origin in the energy method used 
by Diriehlet to address the boundary value problem in classical electrostatics that was 
subsequently named after him. When D is a bounded Lipschitz domain, Bass and 
Hsu [5] constructed the reflecting Brownian motion living on D by showing that the 
so-called Martin-Kuramochi boundary coincides with the Euclidean boundary in this 
case. A general diffusion process on a bounded Lipschitz domain, even allowing locally 
a finite number of Holder cusps, was first constructed by Fukushima and Tomisaki [26]. 
In this work, we use such a Diriehlet form construction in order to derive Feynman-Kac 
representation formulae for the solutions of Neumann, respectively Robin, boundary 
value problems modeling EIT measurements. Probabilistic approaches to both, para¬ 
bolic and elliptic boundary value problems for second order differential operators have 
been studied by many authors, starting with Feynman’s Princeton thesis [20] and the 
article [34] by Kac. The probabilistic approach to the Diriehlet problem for a general 
class of second-order elliptic operators with merely measurable coefficients, even allow¬ 
ing singularities of a certain type, was elaborated by Chen and Zhang [14]; see also 
Zhang’s paper [63]. However, there are only few works that treat Feynman-Kac rep¬ 
resentation formulae for Neumann or Robin type boundary conditions. Moreover, the 
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approaches existing in the literature consider either the Laplacian, see, e.g., [5, 9, 31], 
or non-divergence form operators with smooth coefficients, see, e.g., [22, 49, 10]. For 
the particular case of the conductivity equation on a bounded Lipschitz domain, we 
generalize both, the Feynman-Kac formula for the Robin problem on domains with 
boundary of class for an isotropic C'^’'''-smooth conductivity, 7 > 0, obtained by 
Papanicolaou [49] as well as the representation obtained by Bencherif-Madani and Par- 
doux [10] for the Neumann problem under similar regularity assumptions. While both 
of the aforementioned approaches use stochastic differential equations and Ito calcu¬ 
lus, our approach is based on the theory of symmetric Dirichlet forms, following the 
pioneering work [5] for the reflecting Brownian motion by Bass and Hsu. We derive in 
this work Feynman-Kac formulae for both, the Robin boundary value problem, corre¬ 
sponding to the so-called complete electrode model, as well as the Neumann boundary 
value problem, corresponding to the so-called continuum model. Both formulae are 
valid for possibly anisotropic, uniformly elliptic and uniformly bounded conductivities 
with merely measurable coefficients on bounded Lipschitz domains. During the prepa¬ 
ration of this work we became aware of the paper [15] by Chen and Zhang, where a 
probabilistic approach to some mixed boundary value problems with singular coeffi¬ 
cients is derived. In contrast to our setting, however, the mixed boundary condition 
studied there results from a singular lower-order term of the differential operator. 

Homogenization of reflected stochastic differential equations and partial differential 
equations with Neumann boundary conditions, respectively, in half- space type domains 
have been studied for periodic coefficients in [4, 6 , 59] and for random divergence form 
operators with smooth coefficients in [52]. In contrast to boundary value problems 
with homogeneous Dirichlet boundary conditions, these problems are non-translation 
invariant, which excludes the standard stochastic homogenization approach via the so- 
called environment as viewed from the particle. Employing the Feynman-Kac formula 
in conjunction with a recently obtained invariance principle for reflecting diffusion 
processes associated with random divergence form operators with merely measurable 
coefficients due to Chen, Croydon and Kumagai [13] we provide a homogenization 
result for a stochastic forward problem built on the complete electrode model. Clearly, 
such a result motivates the derivation of stochastic numerical homogenization schemes 
for the approximation of the effective conductivity which are based on simulation of the 
underlying diffusion processes. However, the convergence analysis of such a method 
requires a quantitative convergence result that is stronger than the usual qualitative 
results obtained from the central limit theorem for martingales. As in the case of a 
discrete random walk in random environment, cf. Gloria and Mourrat [27], it turns 
out that the behavior at the bottom of the spectrum of the infinitesimal generator of 
the environment as viewed from the particle process, projected on a suitably chosen 
function, yields bounds on the approximation error. This spectral behavior has been 
the subject of recent interest. Most notably, Gloria, Neukamm and Otto [28] have 
obtained optimal estimates in the discrete case which have been carried over to the 
continuum case by Gloria and Otto [29] . The main difficulty in obtaining such estimates 
for diffusion processes evolving in random media arises from the lack of a Poincare 
inequality for the horizontal derivative in the space of square integrable functions on 
the probability space which corresponds to the random medium. Therefore, in contrast 
to the periodic case, where the Poincare inequality on the torus is available, one can 
not expect a spectral gap in the random case. Still, it has been shown that the bottom 
of the spectrum is sufficiently “thin”. Using these estimates together with a classical 
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argument due to Kipnis and Varadhan [40], we obtain an estimate for the speed of 
convergence of the projected mean-square displacement of the underlying diffusion 
process in a random medium to its limit. Qualitative results of this kind have been 
obtained by Kipnis and Varadhan in the case of discrete random walks in random 
environments and by De Masi, Ferrari, Goldstein and Wick [18] in the continuum 
case, whereas qualitative results in the case of discrete random walks have been proved 
more recently by Gloria and Mourrat [27] and Egloffe, Gloria, Mourrat and Nguyen 
[19]. Finally, we refer to the paper [44] by Mourrat which initiated the idea of using 
the Kipnis and Varadhan argument in order to obtain quantitative results. 

The rest of the paper is structured as follows: We start in Section 2 by briefly intro¬ 
ducing our notation. In Section 3, we recall the modeling of electrode measurements in 
EIT as well as the modeling of random heterogeneous media. Moreover, we introduce 
the stochastic forward problem we are interested in. In Section 4, we describe the con¬ 
struction of reflecting diffusion processes via Dirichlet form theory and in Section 5 we 
derive Skorohod decompositions for two practically relevant classes of conductivities. 
Subsequently, in Section 6, the Feynman-Kac formulae for the deterministic bound¬ 
ary value problems will be derived. Then in Section 7 we study the interconnection 
between Feynman-Kac formulae, stochastic homogenization and stochastic numerics. 
Finally, we conclude with a brief summary of our results. 


2. Notation 

Let D denote a bounded Lipschitz domain in W^, d > 2, with connected complement 
and Lipschitz parameters (ro,co), i.e., for every x € dD we have after rotation and 
translation that dD n B^x^td) is the graph of a Lipschitz function in the first d — 1 
coordinates with Lipschitz constant no larger than cd and D (1 B(x,ro) lies above the 
graph of this function. Moreover, we set := {x G : a; • < 0}, with v = ed the 

outward unit normal on where we identify the boundary of Rl with with 

straightforward abuse of notation. 

For Lipschitz domains, there exists a unique outward unit normal vector v a.e. on 
dD so that the real Lebesgue spaces LP{D) and LP{dD) can be defined in the standard 
manner with the usual norms ll’Hp, p = l,2,oo. The standard inner-products 
are denoted by (•,•) and respectively. The d-dimensional Lebesgue measure is 

denoted by m, the (d — l)-dimensional Lebesgue surface measure is denoted by a and 
I • I denotes the Euclidean norm on R"^. 

By (F, Q, V) we always mean a complete probability space corresponding to a random 
medium. We use the notation ui for an arbitrary element of F and M for the expectation 
with respect to the probability measure V. We use bold letters to denote functions 
on (r,C/,7^), while we use italic letters for the corresponding realizations on R"^ x 
F. The canonical probability space corresponding to diffusion processes evolving in 
a deterministic medium starting in x is denoted P^,) and the expectation with 

respect to P^, is denoted E^,. If the process is evolving in a random medium, we indicate 
this with a superscript uj for the probability measure, i.e., the measure P“ corresponds 
to the particular realization w of the medium. Finally, the product probability space 
corresponding to the annealed measure P:=PPq onfl :=Fxfl, which is obtained 
by integrating with respect to the measure Pq and subsequent averaging over the 
realizations of the random medium, is denoted (fl, P). The expectation with respect 
to P is denoted E. 
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All functions in this work will be real-valued and derivatives are understood in 
distributional sense. We use a diamond subscript to denote subspaces of the standard 
Sobolev spaces containing functions with vanishing mean and interpret integrals over 
dD as dual evaluations with a constant function, if necessary. For example, we will 
frequently use the spaces 

:= e : (0,1)^^ = o} 

and 

HliD) :={<^ei?ip):(</),l)=0}. 

Moreover, we will frequently assume that dD is partitioned into two disjoint parts, 
diD and d 2 D. We denote by Hq{D U diD) the closure of C^{D U diD), the linear 
subspace of C°°{D) consisting of functions (p such that supp(^) is a compact subset of 
D U diD, in H^\d). Moreover, we define the Bochner space A^(r; Hq{D U diD)) 

= {0 : r ^ U d^D) : dT’H < 00 }, 

see, e.g., [3] for properties of this space. 

For the reason of notational compactness, we use the Iverson brackets: Let S' be a 
mathematical statement, then 




1, if S is true 

0 , otherwise. 


We also use the Iverson brackets [x G B] to denote the indicator function of a set B, 
which we abbreviate by [B] if there is no danger of confusion. 

In what follows, all unimportant constants are denoted c, sometimes with additional 
subscripts, and they may vary from line to line. 


3. Electrical impedance tomography forward problems 

3.1. Modeling of electrode measurements. We assume that the, possibly anisotropic, 
conductivity is defined by a symmetric, matrix-valued function «:£)—>■ with 

components in L°°{D) such that k is uniformly bounded and uniformly elliptic, i.e., 
there exists some constant c > 0 such that 

(2) ^ • k{x)^ < c|^p, for every ^ and a.e. x G D. 

The forward problem of electrical impedance tomography can be described by differ¬ 
ent measurement models. In the so-called continuum model, the conductivity equation 
( 1 ) is equipped with a co-normal boundary condition 

(3) QkvU := Kv ■ Vitjau = / on dD, 

where / is a measurable function modeling the signed density of the outgoing current. 
The boundary value problem (1), (3) has a solution if and only if 

(4) (/) l)a_D = 

Physically speaking, this means that the current must be conserved. Given an appro¬ 
priate function /, the solution to (1), (3) is unique up to an additive constant, which 

_ /2 

physically corresponds to the choice of the ground level of the potential. If / G iL* , 
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then there exists a unique equivalence class of functions u € H^{D)/M. that satisfies 
the weak formulation of the boundary value problem 


kVu • Vw da: = (/, v\dD) 


ID 


dD 


for all V € H^{D)/R, 


where v\qd '■= "fv and 7 : H^{D)/M. — >■ = (iJo is the standard 

trace operator. Note that we occasionally write v instead of v\dD for the sake of 
readability. 

In practical EIT measurement, a number of electrodes, denoted Ei,...,E!^ C dD^ 
are attached on the boundary of the object D. These electrodes are modeled by dis¬ 
joint surface patches given by simply connected subsets of dD, each having a Lipschitz 
boundary curve. The most accurate forward model for real-life EIT is the so-called 
complete electrode model which takes into account the fact that during electrode mea¬ 
surements there is a contact impedance caused by a thin, highly resistive layer at the 
electrode object interface. It was demonstrated experimentally that the complete elec¬ 
trode model can correctly predict measurements up to instrument precision, cf. [57]. 
For a given voltage pattern U G the boundary conditions for the complete electrode 
model are given by 


(5) -yuloD + gu\dD = f ondD, 

where the functions /, g : dD —>■ R are defined by 

N N 

^ ' 1=1 ^ ^ 1=1 

and the contact impedance z : dD —>■ R is assumed to be a piecewise continuous 
function, with interfaces that are of zero surface measure, satisfying 


0 < Co < z < Cl a.e. on dD. 

For a given voltage pattern U G satisfying the grounding condition 

N 

(7) 

1=1 

the equations (1) and (5) define the electric potential u G H^{D) uniquely, cf. [57] and 
the variational form of the boundary value problem (1), (5) reads as follows: Given 
U G R^ satisfying (7), find u G H^{D) such that 

( 8 ) / KVu-Vvdx + {gu\dD,v\aD)gD = {f,v\oD)9D ioT all V G H^{D). 

J D 

Knowledge of u yields the corresponding electrode current vector J G R^ via 

( 9 ) Ji = [ df,^uda(x), 1<1<N. 

J El 


3.2. A stochastic forward problem. The basic geometric setting of the stochastic 
problem we are interested in is as follows: Assume that the model domain is given by 
the lower hemisphere 


£> := B(0,i?)nRl, R>Q 


and that dD is partitioned into two disjoint parts, namely the accessible boundary 
diD := dDr\R‘^~^ and the inaccessible boundary d 2 D := dD\diD, respectively. Such 




From Feynman-Kac formulae to stochastic homogenization in EIT 


7 


a setting is found for instance in geophysical applications, where measurements can 
only be taken on the surface, cf., e.g., [55]. 

Let (r, Q, V) be a probability space and let © : L —>■ L denote an ergodic d- 
dimensional dynamical system, i.e., a family of automorphisms {©x, a: € which 
satisfies the following conditions: 

(i) The family {©x,a: G R'^} is a group, i.e., ©o = id and 

©x+y = ©x®y for all x,y G 

(ii) the mappings ©x : T —T, x G preserve the measure P on F, i.e., for every 
B G G, ®xB is 7^-measurable and 

V{®xB)=V{B)- 

(hi) for every measurable function cf) on {F,G,'P), the function (x,w) cf>{&xU!) is 

a measurable function on (R"^ x T, B{W^) ^G^cuxT), where 0 G denotes 

the sigma-algebra generated by the measurable rectangles; 

(iv) the family {©x,a: G R'^} is ergodic, i.e., 0(©xw) = cf>{u}) for all x G R"^ and 
P-a.e. O' G r implies cf> = const V-a.e. 

Throughout this work we assume that the conductivity random field 

{k{x,lo), {x,u}) G R'^ X T} 

is the stationary extension on R'^ x Tof some function k G L^(r;that is, 

(10) ix,uj) I—>■ k{x,uj) = k(©xw). 

Note that if k can be written in the form (10) with a dynamical system {©x,x G R'^} 
which satisfies (i)-(iii), then it is automatically stationary with respect to V, i.e., for 
every finite collection of points x*^®\ * = 1, k, and any h G R'^ the joint distribution 


/t(x^^^ -1- h,ijS),. 

under V is the same as that of 

.., /t(x^^^ -1- h, w) 

/t(x(^\w),.. 

.,k{x^^\u)). 


Even if it is not explicitly stated, we always assume that the conductivity random held 
may be written in the form (10), where the underlying dynamical system {©x,x G 
R'^} satishes conditions (i)-(iv). Moreover, we will explicitly state if the conductivity 
random Held {k{x,uj), {x,uj) G R'^ x T} satishes one of the following assumptions: 

(Al) K G L^(r; L°°(R‘^; R‘^^‘^)) and the random held is strictly positive and uni¬ 
formly bounded, that is, there exists a constant c > 0 such that for every 
^ G R'^ and a.e. x G R'^ 

r^UJ G r : < C • k{x,uj)^ < cj^pj = 1. 

(A2) {K(x,a;), (x,a;) G R*^ x T} satishes the spectral gap property, cf. [29]: There 
exist constants p > 0 and r < oo such that for all measurable functions on 
{k : R'^ ^ {/to G : J/to^] < lCl,cl^l2 < ^ •/toC for all ^ G R'*}} 

YcI)<-m[ ( osc«|3 cf>) dx, 

P JkP ^ ' ' 


Petteri Piiroinen and Martin Simon 


where we have set 

(oSC^lBC^.r) :=SUp|0(k) : k € Q, k\Kd\B{x,r) = 

- inf j </>(k) -.ken, = K 


}■ 


To account for the highly heterogeneous properties of the background medium, the 
latter is modeled using the conductivity random field with appropriate scaling by a 
small parameter £ > 0, i.e., 

Ke{-, •) : X r — Ke(x, uj) := K(a;/e, w). 

If the correlation length of the conductivity random field k is, say 1, then the correlation 
length of the scaled version is of order e and for e <C 1 we obtain thus a rapidly 
oscillating random field. 

3.3. Stochastic forward and inverse problem. Now let us introduce a stochastic 
forward model based on the complete electrode model: We search for a random held 
{u^(x,uj), (x,uj) G D X r} for the electrical potential with G L'^(T-, Hq{D U diD)) 
such that the stochastic conductivity equation 

(11) V • {K^\7ug) =0 in I? X r 

subject to the boundary conditions 

-VusldiD + gUeldiD = f ondiDxT 
UE\d 2 D = 0 on d 2 D X r 

is satished V-sl.s. The variational formulation of the forward problem is to hnd G 
L'^{T; Hq{D U diD)) such that 


( 12 ) 


(13) 


I J KeVUe • Vu dx + {gUe, v )} = ^ (/, v) 


for all V G T^(r; Hl{D\JdiD)). For a given voltage pattern U G R^ the corresponding 
measurement data is given by the random current pattern J{£,uj) = (Ji(e,a;),..., Jn (s) w))^, 
dehned for P-a.e. a; G T by 


(14) 


■Wue{x,uj)\oiDda{x), 1 = 1, 

\^l\ JEi 


Due to the assumption (Al), the well-posedness of the variational formulation follows 
from a straightforward application of the Lax-Milgram theorem. Moreover, standard 
arguments from measure theory show that the solution to the stochastic forward prob¬ 
lem (11), (12) also solves (13), cf. [3]. 


4. Reflecting diffusion processes 

In his seminal paper [23], Fukushima established a one-to-one correspondence be¬ 
tween regular symmetric Dirichlet forms and symmetric Hunt processes, which is the 
foundation for the construction of stochastic processes via Dirichlet form techniques. 
Therefore we assume that the reader is familiar with the theory of symmetric Dirichlet 
forms, as elaborated for instance in the monograph [25]. 

Let us consider the following symmetric bilinear forms on L^(D): 

(15) £{v,w) := f k\7v{x) ■ \7w{x) dx, v,w £ !){£) := H^{D) 

JD 
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and for the particular case k = 1/2, which is of special importance, we set 


(16) £^^{v,w):=l- f Vv{x)-Vw{x)dx, v,w € V{£^^) := H\D). 

^ .1 n 


The pair {£,'D{£)) defined by (15) is a strongly local, regular symmetric Dirichlet 
form on L?{D). In particular, there exist an ^-exceptional set Af C D and a conserva¬ 
tive diffusion process X = {£1, X, {Xt, t > 0}, P^,), starting from x G D\J\f such that X 
is associated with (f, 'D{£)). Without loss of generality let us assume that X is defined 


on the canonical sample space ft = (^([O, oo); D). It is well-known that the symmetric 


Hunt process associated with (16) is the reflecting Brownian motion. Therefore, we 
call the symmetric Hunt process associated with (15) a reflecting diffusion process. 

Let us briefly recall the concept of the boundary local time of reflecting diffusion 
processes, see, e.g., [31, 49, 10], which will be crucial for the subsequent derivation 
of the Feynman-Kac formulae. If the diffusion process is the solution to a stochastic 
differential equation, say the reflecting Brownian motion, then the boundary local 
time is given by the one-dimensional process L in the Skorohod decomposition, which 
prevents the sample paths from leaving D, i.e., 



(17) 


P2,-a.s. for q.e. x € D. This boundary local time is a continuous non-decreasing 
process which increases only when Xt G dD, namely for alH > 0 and q.e. x G D 



Pa,-a.s. and 



Although the reflecting diffusion process associated with (15) does in general not admit 
a Skorohod decomposition of the form (17), we may still define a continuous one¬ 
dimensional process with these properties. More precisely, by the Lipschitz property 
of dD, we have that Dr\B{x,rD) = {{x,Xd) : Xd > 'y{x)}r\B{x,rD) and the Lipschitz 
function 7 is differentiable a.e. with a bounded gradient. In particular, we have for 
every Borel set B C dD fl B{x, riff) that 

cr(H) = [ (l + |V7(i)p) ^ di 

and a straightforward computation yields that the Lebesgue surface measure cr is a 
smooth measure with respect to {£,D(£)) having finite energy, i.e.. 



where we have used the inner product £i{-, •) := £{■, ■) + (•, •). 

Definition 4.1. The positive continuous additive functional of X whose Revuz mea¬ 
sure is given by the Lebesgue surface measure a on dD, i.e., the unique L G such 
that 


(18) 
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for all non-negative Borel functions (j) and all a-excessive functions is called the 
boundary local time of the reflecting diffusion process X. 


The rest of this section is devoted to showing that the ^-exceptional set Af is actually 
empty. Therefore, we consider the non-positive definite self-adjoint operator {C,'D{C)) 
associated with the Dirichlet form {£,'D{£)). That is, for v G 'D{C) we have 

(19) w) = £{v, w) for all w € £>{£) 


and the domain of C is given by 

£>{£) = G TA{£) '■ 3^ G L^{D) s.t. £(v,w) = J cjiwdx Ww G T>{£)^. 


In order to refine the reflecting diffusion process X to start from every x G D, we exploit 
the connection between the strongly continuous sub-Markovian contraction semigroup 
{Tt,t > 0} on Lf{D) and the evolution system corresponding to (£,T’(£)), see, e.g., 
the monograph [50]. Namely, for every vq G Lf{D)^ the trajectory v : (0,r) — >■ H^{D), 
v{f) = Ttvo belongs to the function space 

{</> G L\{Q,T)-H\D)) : 0 G L\{Q,T)-H-\D))} 


and is the unique mild solution to the parabolic abstract Cauchy problem 


( 20 ) 


V + Cv = 0 in (0, T) 
u(0) = Uq. 


This is equivalent to the variational formulation 

(21) — f {v{t),w) (p{t) dt + [ {Cv{f),w) ipft) dt — {vq,w) (p{0) = d 

do do 

for all w G H^{D) and all ip G C'“([0,r)). Moreover, Tt is known to be a bounded 
operator from L^{D) to L°°{D) for every t > 0. Therefore, by the Dunford-Pettis 
theorem, it can be represented as an integral operator for every t > 0, 

(22) Tt(j){x) = [ p{t,x,y)(j){y)dy for every </> G L^iD), 

J D 

where for alH > 0 we have pit,-,-) G L°°{D x D) and p{t, -,-) >0 a.e. We call the 
function p the transition kernel densityoi X. 

The following proposition adapts a well-known result for diffusion processes on K.'^, 
cf. [58], which follows from the famous De Giorgi-Nash-Moser theorem, to the case of 
reflecting diffusion processes on D. The key idea of the proof is the following extension 
by reflection technique from [62, Section 2.4.3]: We extend the solution to a parabolic 
problem by reflection at the boundary. Then we show that this extension again solves a 
parabolic problem so that we can apply the interior regularity result due to De Giorgi, 
Nash and Moser. See also the article [46] by Nittka, where such a technique is applied 
to elliptic boundary value problems. 

Proposition 1. p G C°’^((0,T] x D x D) for some 6 G (0,1), i.e., for each fixed 
0 < to < T, there exists a positive constant c such that 

(23) \p{t 2 ,X 2 ,y 2 ) -p{ti,xi,yi)\ < c{y/t 2 - h -b \x 2 - a;i| -b \y 2 - yi\)^ 

for all to < ti < t 2 < T and all (xi,yi), {x 2 ,y 2 ) G D x D. Moreover, the mapping 
1 1—>■ pit, - ,- ) is analytic from (0, oo) to x D). 
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Proof. First note that Nash’s inequality holds for the underlying Dirichlet form {£, {D)), 

i.e., there exists a constant ci > 0 such that 

||^||2+4/d ^ ci{£{v, v) + ||u|||)||ti||f^‘^ for all v € H^{D). 

This is a direct consequence of the uniform ellipticity (2) and [5, Corollary 2.2], where 
Nash’s inequality is shown to hold for the Dirichlet form {£^^,H^{D)) for a bounded 
Lipschitz domain D. Analogously to the proof of [5, Theorem 3.1], it follows thus 
from [12, Theorem 3.25] that the transition kernel density satisfies an Aronson type 
Gaussian upper bound 

(24) p(t,a;,2/) < cif'^/^exp 

for all t < 1 and all {x,y) G D x D. In particular, supQ.^j<i||p(t, •, •)||oo is finite 
and hence by the interior Holder continuity obtained from the De Giorgi-Nash-Moser 
theorem, cf. [45, 58], the estimate (23) is true for all {x2,y2) satisfying 

d{xi,dD),d{yi,dD) > ca, * = 1,2, for some constant ca > 0 and all to < ti < t2 < 1. 
Note that by the semigroup property the Chapman-Kolmogorov equation holds, i.e., 

(25) p(ti+t2,x,y) = / p{ti,x,z)p{t2,z,y)dz 

Jd 

for every pair ti,t2 > 0 and a.e. x,y G D. In particular, for fixed y G D the function 
V := p{-,-,y) is the unique solution to (20) with initial value vo ■= p{0,-,y) G L^{D). 
Now let 2 G dD so that by the Lipschitz property of dD we have after translation and 
rotation B{z, ro) fl D = {(i, Xd) G B(z, rn) '■ Xd > 7(i)} and B{z, ro) fl dD = {a; g 
B{z,ro) : Xd = 7(5:)}, where we have introduced the notation x = (cci, ..., Xd-i)^- 
Let us furthermore introduce the one-to-one transformation \I'(a;) := {x,Xd — 7(5)) 
which straightens the boundary B(z, ro) C\dD. d' is a bi-Lipschitz transformation and 
the Jacobians of both and are bounded with bounds that depend only on the 
Lipschitz constant cd- Since v is the solution to (20) with appropriate initial condition, 
the function v := v{-, 4'“^(-)) must satisfy the variational formulation of the following 
parabolic problem in D{z,rD) '■= d>{B{z,rD) H D) 




v(t)w dx dt 


— > / ip{t) / kijdiv{t)djW dx dt 

^ Jb{z,rD) 


i,j=l 


— 1 ^( 0 ) / vowdx 

J b{z,ro) 


for all w G C^{D{z,rD)) and all (p G C)]’°([0,T)). The coefficient k is obtained via 
change of variables and it is bounded and uniformly elliptic by the boundedness of the 
Jacobians of 4' and 4'“^, respectively. Now we use reflection at the hyperplane {{y, 0)} 
via the mapping p{x) := (5, —Xd) which yields that the function v{-,p{-)) satisfies 
the variational formulation of a parabolic problem on p{D(z.,rr))). Summing up both 
variational formulations on D{z,r£i) and on p{D{z,r£i)), respectively, we obtain that 
the function 

2.') xGD{z,rn) 

\v{t,p{x)), xGp{D{z,rD)) 

satisfies the variational formulation of a parabolic problem in D(z,ro) U p(D{z,ru). 
By the interior Holder estimate for v, together with the fact that we may choose 
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C3 = td/^cd, we obtain thus 

\p{t2,X2,^~^{y2)) -p{ti,xi,^~^{yi))\ < ci(Vi2 - ti + \y2 - yi|)“" 

for all to < *1 < ^2 < 1 and yi,y 2 G {{x,Xd) : |i| < C3, Xa G (0,r£)/4)}. As 4* 
is bi-Lipschitz, for fixed x, the mapping (t,2/) 1—>■ p{t,x,y) is Holder continuous in 
(to,l] X (H(z, C3) n D) and by symmetry of the transition kernel density the same 
holds true for the mapping (t, x) 1—>■ p{t, x, y) for fixed y. Finally, the hrst assertion on 
(to,!] X D X D follows due to compactness of dD and its generalization to arbitrary 
T > 0 is obtained after repeatedly applying the Chapman-Kolmogorov equation. 

The second assertion follows by the fact that the semigroup {Tt,t > 0} extrapolates 
to a holomorphic semigroup on L‘^{D). More precisely, the semigroup possesses a 
holomorphic extension to the sector Eg := {re“ : r > 0, |a| < 9} for some 9 G (0, ^], 
cf., e.g., [50]. Let 0 < to <T and set 

2^) '■= {z G C : z — to G Eg, jzj < T}. 

By the Holder continuity of p, the set {p{z,-,-) : z G T,g{to,T)} is a bounded subset 
of C^’^{D X D). Moreover, the family of functionals obtained form integration against 
the functions [Bi\{x)[B2\{y) for measurable Bi, B2 C D form a separating subspace of 
{C°’^(D,D)y, i.e., for k G C°’^(D x D) 

/ k{x,y)[Bi]{x)[B 2 ]{y) dxdy = 0 for all measurable Hi,i?2 C D 

Jdxd 

implies that k = 0. As the mapping 

z^ {T^[Bi],[B 2 ]) = f p{z,x,y)[Bi]{y)[B 2 ]ix)dxdy 
JDxD 

is holomorphic for all z G Eg, the mapping z 1—>■ p{z, •, •) is holomorphic from Eg(to, T) 
to C^’^{D X D) by [1, Theorem 3.1]. Since to and T were arbitrary, the assertion is 
proved. □ 

By [24, Theorem 2], the existence of a Holder continuous transition kernel density 
ensures that we may refine the process X to start from every a; G H by identifying the 
strongly continuous semigroup {Tt.,t > 0} with the transition semigroup {Pt,t > 0}. 
In particular, if v is continuous and locally in H^{D), the Fukushima decomposition 
holds for every x € D, i.e., 

(26) v{Xt) = v(Xo) + My + Xy, for all t > 0, 

Pa:-a.s., where M'" is a martingale additive functional of X having finite energy and 
N'’ is a continuous additive functional of X having zero energy. 

Moreover, both M'" and N'" can be taken to be additive functionals of X in the 
strict sense, cf. [25, Theorem 5.2.5]. 

Finally, note that the 1-potential of the Lebesgue surface measure a of dD is the 
solution to an elliptic boundary value problem on a Lipschitz domain with bounded 
data. By elliptic regularity theory, cf., e.g., [30], this solution is continuous, implying 
that the boundary local time L exists as a positive continuous additive functional in 
the strict sense, cf. [25, Theorem 5.1.6]. 
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5. Skorohod decompositions 


In this section, we derive Skorohod decompositions of the reflecting diffusion process 
X for two practically relevant special cases, namely local Lipschitz conductivities and 
isotropic piecewise constant conductivities. 

The assertion of the following proposition is already covered by [26, Theorem 2.3]; 
we include a proof for the sake of self-containedness. 

Proposition 2. Let k G be a symmetric, uniformly bounded and uni¬ 

formly elliptic conductivity. Then the reflecting diffusion process X admits the follow¬ 
ing Skorohod decomposition 

VK{Xs)ds- [ K{XsHXs)dLs, 

Jo 

¥x-a.s., where B : D ^ denotes the positive definite diffusion matrix satisfying 

= 2k, W is a standard d-dimensional Brownian motion and L is the boundary local 
time of X. 


(27) 


Xt=x- 


BiXs)dWs 


Proof. We have shown in Section 2, that the Fukushima decomposition holds with a 
unique martingale additive functional AT" in the strict sense and a unique continuous 
additive functional N" in the strict sense. Let us first compute the energy measure of 
M". For v,w G T){6) we obtain 


'D 


w{x)dp.(^M-){x) = 


lim i [{v{Xt) - v{xo)f'] w{x) dx 

^0+ t Jjj ^ 

= lim - j {Ttv'^{x) — 2v{x)Ttv{x)-\-v'^{x))w{x) dx 
*->0+ t J 

= lim - [ v(x)w{x)(v(x) — Ttv(x)) dx 
t^o+ t Ju 

— lim - j v'^{x){w{x)—Ttw{x))dx 

t->0+ t Jjj 

= 2£{vw, v) — £(v^,w) 

= 2 kVv{x) ■ Vv{x) w{x) dx, 

Jd 


which yields the energy measure 


d 

dfJ.{M-){x) = 2 E Kij{x)div{x)djv{x) dx 

so that the predictable quadratic variation of M" is given by 

i-t d 

(28) {xr)^ = 2 ^/c,,(x,)a,u(x,)aju(x,)ds. 

JO 1 

i,j=i 


Using the coordinate mappings 4>i{x) := Xi, i = l,...,d, on D yields that AT’^ is a 
continuous martingale additive functional in the strict sense with covariation 


{AI^fAT^^)^ 


= 2 / Kij{Xs)ds, 

Jo 
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Pa^-a.s. A standard characterization of continuous martingales, cf., e.g., [?], yields that 

(29) Mj" = f iB{X,)Vv{Xs)) ■ diy,, 

Jo 

Px-a.s., where B : D ^ denotes the positive definite diffusion matrix satisfying 

B^ = 2k and W is a standard d-dimensional Brownian motion. 

Now let us consider the continuous additive functional N'^. Again using the coordi¬ 
nate mappings on D, we obtain from Green’s formula that 


d „ 

= 2, / Kij(x)dj'w{x) dx 

d ^ d ^ 

= — / djKij{x)w{x) dx + 2, / Kij{x)yj{x)w{x) dcj{x) 

,= 1 Jd JdD 


for all V G H^{D) fl C{D). That is, is associated with the signed Radon measure 

d d 

— djKij da; -f Kij{x)vj(x) da(x) 

i=i i=i 

and by the fact that the unique positive continuous additive functionals in the strict 
sense having as Revuz measure the Lebesgue measure and the Lebesgue surface mea¬ 
sure, respectively, are given by the constant additive functional t and the boundary 
local time Lt, we have shown that for every x € D 

d pt d pi 

(30) Np=y^ a,Ky(X,)ds-V / K,jiXs)iyjiXs)dLs, 

Pa,-a.s. Substitution of (29) and (30) in the Fukushima decomposition for the coordi¬ 
nate mappings finally yields the Skorohod decomposition (32). □ 


Now let us turn to the case of isotropic piecewise constant conductivities and for 
simplicity of the presentation let us consider a simplistic two-phase medium, where 


(31) 


Ki, X € Di 
K2, X G D 2 , 


with constants ki,K 2 > 0 and D is a simply connected bounded Lipschitz domain 
which consists of two disjoint subdomains such that Di = D\D 2 . We assume that D 2 
is a simply connected Lipschitz domain, v is the outer unit normal vector on dD and 
the outer unit normal vector on dD 2 with respect to I?2- 


Definition 5.1. The positive continuous additive functional of X whose Revuz 
measure is given by the scaled Lebesgue surface measure (ki -I- K 2 )ct on dD 2 is called 
the symmetric local time of the reflecting diffusion process X at dD 2 . 


Remark 1. The term “symmetric” comes from the fact that in the one-dimensional 
case is the local time defined by the Tanaka formula with the convention sign(O) = 0, 
which is called the symmetric local time, see [?]. In this case we have 
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Proposition 3. Let k he given by (31 ). Then the reflecting dijfusion process X admits 
the following Skorohod decomposition 

(32) Xt=x + [ [ ,y{Xs) dL° [ iy{Xs) dLs, 

Jo 1^1 + «2 Jo Jo 

Tx-a.s., where W is a standard d-dimensional Brownian motion, LP is the symmetric 
local time of X at dD 2 and L is the boundary local time. 


Proof. Repeating the computations from the proof of Proposition 2 yields for the mar¬ 
tingale additive functional the predictable quadratic variation 


= 2 



K(X,)|Vz;(X,)|2ds, 


implying that 


Pa;-a.s. for every x G D, 


Mf = [ y'2K{X,)Xv{X,)) ■ dlP„ 


c-a.s. for every x € D, 


Jo 

where W is a, standard d-dimensional Brownian motion. By Green’s formula we have 
for all w G PiS) n C{D) 


£{v,w) = Ki Vn • Vrc dx-I-K2 / Vn-Vu>da; 

J J D 2 

= — vkAv dx — {ki — K2) / d,yvwda{x) 

Jd JdD2 

+Ki / dvVW d< 7 {x). 

JdD 


Using the coordinate mappings on D we obtain that is associated with the signed 
Radon measure 

-(ki - K2)v{x) daldD^ix) + Kiv{x) da\QD{x). 

The assertion follows by the fact that the unique positive continuous additive function¬ 
als in the strict sense having as Revuz measure the (scaled) Lebesgue surface measure 
on dD 2 and dD, respectively, are given by the symmetric local time at dD 2 and 
the boundary local time Lt. □ 


6. Feynman-Kac formulae 

In this section, we derive the Feynman-Kac formulae for both the continuum model 
and the complete electrode model. Afterwards, we will obtain, as a corollary, a 
Feynman-Kac formula for the mixed boundary value problem corresponding to the 
stochastic anomaly detection problem introduced in the first chapter. Compared to 
the earlier works [31, 49, 10] on Feynman-Kac formulae, the main difficulty in deriving 
these formulae in our particular setting comes from the lack of Ito’s formula for general 
reflecting diffusion processes. 

The rest of this subsection is devoted to providing some auxiliary lemmata. 

Lemma 6.1. The transition kernel density p approaches the stationary distribution 
uniformly and exponentially fast, that is, there exist positive constants ci and C 2 such 
that for all {x,y) € D x D and every t > 0, 

(33) \p{t,x,y) - \D\~'^\ < ciexp(-C2t). 
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Proof. Fix X € D. Then by the Chapman-Kolmogorov equation 

p{t,x,x) - \D\~^ = [ p{t/2,x,y)p{t/2,y,x)dy -\D\~'^ 

J D 

= [ {p{t/‘^:X,y))^dy-\D\-^ 

J D 

= [ (P(V2, X, y) - \D\-^f dy > 0, 

J D 

where we have used that Jjjp{t,x,y) dy = 1 for alH > 0 and x G D. Moreover, we 
have by the analyticity of the mapping 1 1—^ p(t, x, x), cf. Theorem 1, and the fact that 
p solves a parabolic boundary value problem with homogeneous Neumann boundary 
condition 

d d /* 

— {p{t,x,x) - = -^ J pit/‘^,x,y)p{t/2,y,x)dy 


< 

< 


^ J^ipit/‘^,x,y)fdy 
- / i^{y)"^vP{t/2,x,y) ■Vyp{t/2,x,y)dy 

JD 

/ '^yPitl‘2,x,y) ■Vyp{t/2,x,y)dy 
JD 

{ [ {p{t/2,x,y)rdy-\D\-A 


m 

-1 


—c c 


D 


= -c ^Cjj^{p{t,x,x) - \D\ 
where we have used (2) and the Poincare inequality 


(f) —\D\ ^ f (j){x) dx 
JD 


< C£)||V(^||2 ioT all (j) G {D). 


Integration of the inequality from above yields a diagonal estimate, i.e., there exist 
positive constants ci and C2 such that 

0 < p{t, X, x) — \D\~^ < Cl exp(—C2t) for every t > 0. 


Now, the assertion follows from the Cauchy-Schwarz inequality and the fact that, by 
the computations from above, we may write the expression \p{t^x,y) — in the 

form 

[ {p{t/2, X, z) - \D\-^){p{t/2, z, y) - \D\-'^) dz 

JD 

< [J {p{t/‘^,x,z) - dz'^ ^ (^J {p{t/2,y,z)-\D\-'^fdz'j 

< (^p{t/2,x,x) ' [p(tl2,y,y)-\D\~^'^ 

□ 


Lemma 6.2. Let Kij G 1 < i,j < d. Then the set 

(34) V^{D) := {(j) : (j) G = 0 a.e. on dD} n H\D) 

is dense in H^{D). 
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Proof. Diagonalizing the operator {C,'D{L)) corresponding to the conductivity k, we 
obtain an orthonormal basis k gN} of L‘^{D) and an increasing sequence (Afc)fcgN 
of real positive numbers which tend to infinity such that for every fc € N, (/)fc € H^{D) is 
a weak solution of the corresponding eigenvalue problem with homogeneous Neumann 
boundary condition. Note that the inner product £’i(-, •) is equivalent to the standard 
inner product on H^{D). Note that Vn{D) contains the linear span of {(j)k,k € N} 
by elliptic regularity so it is enough to show that the linear span of eigenfunctions is 
dense. Therefore, let V' S H^{D) such that £i{(j)k,'f’) = 0 for every fc S N, then 

0 = / nVcfk ■'^'ipdx + {(fk.i’) = i>^k+ f) {4>k,'ip) ■ 

Jd 

Hence it follows = 0 for every k G N and the fact that {4>k,k G N} is an 

orthonormal basis of Lf{D) implies ip = 0 which proves the assertion. □ 


Lemma 6.3. For every x G D and every bounded Borel function cp on dD we have 


(35) 


E, 


: f 4>{Xs)dLs= [ [ p{s,x,y)cp{y)da{y)ds forallt>0. 

Jo Jo JdD 

Proof. First, the expression (35) is well-defined since the boundary local time of X 
exists as a positive continuous additive functional in the strict sense. Without loss of 
generality we may assume that (p is non-negative. It follows from [25, Theorem 5.1.3] 
that the Revuz correspondence (18) is equivalent to 

ip{x)E^ f p{Xs)dLsdx = f f p{y)Tsip{y)d(j{y)ds 
Jo Jo JdD 

= / / / p{s,y,x)p{y)da{y)dsdx 

Jd Jo Jod 

for every t > 0 and all non-negative Borel functions ip and p, where we have used 
Fubini’s theorem in the second line. As this holds for every non-negative Borel function 
Ip, we deduce 


D 




ft rt 

P{Xs)dLs = 


p{.s,x,y)p{y) dcr{y) ds a.e. in D. 

10 Jo JdD 

To obtain the assertion everywhere in D, fix an arbitrary xo G D and consider for 
to > 0 the integral 

ft-to 


/ PiXs)dL, 

Jto 


j p(to,a;o,a;)E2,| J p{Xs)dLs'^dx 

I p{to,xo,x)( I I P{s,x,y)p{y)da{y)ds) dx 
Jd ^Jo JdD ' 


n p{s,xo,y)p{z) da{y) ds, 

'D 


where we have used the Markov property of X. Now let {tk)kGH denote a positive 
sequence which monotonically decreases to zero as fc —>■ oo. By the computation from 
above we have for every x G D 


E,r. I P{Xs)dLs = 

'0 


p{s,x,y)p{y)da{y)ds+E,^ p{Xs)dLs 


tk JdD 
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The claim follows by the facts that 4> is bounded and ^xLt^ goes to zero as fc —)► oo 
which follows from monotonicity and continuity of the local time and the property 
Lq = 0 P^^-a.s. for every x G D. □ 


6.1. Continuum model. The main result for the continuum model (1), (3) is the 
following theorem. 


Theorem 6.4. Let f be a bounded Borel function satisfying (/, = 0. Then there 

is a unique weak solution u G C{D) fl Hf^D) to the boundary value problem (1), (3). 
This solution admits the Feynman-Kac representation 

(36) u{x) = lim Ea; f f{Xs)dLs for all x G D. 

t-)-oo Jq 


Proof. The existence of a unique normalized weak solution u to (1), (3) is guaranteed 
by the standard theory of linear elliptic boundary value problems. Let us set 

Ut{x):=Ex / f{Xs)dLs and Uao{x) ■= lim ut(a;), x G D, 

Jo 

respectively. From the occupation formula (35) and the compatibility condition (4), it 
follows immediately that 

Mx) = [ f {pis,x,y) - \D\-^)f{y)da{y)ds for all a; G D. 

Jo JdD 


By Lemma 6.1 the convergence towards the stationary distribution is uniform over D, 
in particular, 


(37) Uoo{x)= / {p{s,x,y) -\D\ ^)f{y)d<7{y)ds for all x G T>. 

Jo JdD 

It follows from (37) together with the Holder continuity shown in Proposition 1 and 
the Aronson type upper bound (24) that Uoo is in C{D). Moreover, by the facts that 
p is the transition kernel density of a reflecting diffusion process and / is bounded, 
Fubini’s theorem yields 


/ ( / / P{s,x,y)f{y)d(j{y)ds)dx 

Jd ^Jo JdD ' 




/(2/)dcr(y) ds = 0, 


i.e., Uoo G Hl{D). 

Now, let us use the following regularization technique in order to show u = Uoo- Let 

(K^^^)fcgN denote a sequence of smooth conductivities with components in C°°{D) such 
(k) 

that for 1 < i, j < d, nld Kij a.e. as fe —?► oo. Let us consider the Dirichlet forms 
(fW, 771(11)), fc G N, with 


:= f k’^^Vv-Vwdx 
Jd 


and the associated reflecting diffusion processes . By Proposition 2, we obtain the 
Skorohod decomposition 

X^’"'>=x+[ a('=)(xW)ds+ / B^'^\Xs)dWs- [ dL^^\ 

Jo Jo Jo 

where VF is a standard d-dimensional Brownian motion, a] ^ := i = 

1, and the matrix satisfies Let us define in the same 
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manner as Ut and u^^\x) := limt_>oo x € D. We show that is the unique 

weak solution of the elliptic boundary value problem 

fv-(KWVuW) = 0 inD 

= f ondD 

in the Sobolev space Hl{D). For test functions v £ Vi^(k){D), we may apply Ito’s 
formula for semimartingales to obtain 

)) = vix) + E, / V • ds. 

Jo 

By Fubini’s theorem, this is equivalent to 

tI;'^'^v{x) — v{x) = f f p^^\s,x,y)W ■ {K^'^^Vv{y)) dy ds, 

Jo Jd 

where we have used the superscript “(fc)” for the semigroup and transition kernel 
density, respectively, corresponding to Multiplication with /, integration over 

dD and another change of the orders of integration finally yield 

/ - v{y)) dcr(y) = (wf \ V • . 

JdD ' ' 

Since —>■ vfJJ and \D\~^ jj^vdx, both uniformly on D, as t ^ oo, we 

have 

(zrW,V-(«;WVz;))=-(/,u)9^, 

where we have used the expression (37) with instead of p for As this holds 
true for every v € V),(ic)(D), must be the unique normalized weak solution to the 
boundary value problem by a density argument. 

Next, we show the convergence of the sequence as fc —>• oo towards 

u € Hl{D), the unique solution to (1), (3). From our assumptions on the sequence 
it is clear that —>■ nnj, ^ < i,j < d, in L‘^{D) as fc —>■ 0, which im¬ 
plies G-convergence of the sequence of elliptic operators on D towards £, cf. 

[64]. That is, for any (j) G {Hl{D))', the solutions G Hoi^) ^^e homogeneous 
Dirichlet problem 

(38) V • = (/) inD 

satisfy ^ re in Hq{D) as fe — >■ oo and ^ kVw in L^(D;R'^) as fc —>■ oo, 

where w G Hq{D) is the solution to the homogeneous Dirichlet problem 

V • (kVw) = (j) in D. 

Consider the variational form of the Neumann problem for G i.e.. 


[ ■Vvdx = {f,v)gjj for all u G i7i(D). 

Jd 


As G HI{D) for all fc G N, we have by the Poincare inequality 


||w«||<Cl||Vu('=)||2<C2||/||2. 
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That is, the sequence {u^^'>)kGN is bounded in Hl{D) and we may extract a weakly 
convergent subsequence, still denoted for convenience. We have ^ m in 

HI { D ), as fc —>■ oo, so that it remains to show the convergence of the flows 

^ kVu in L^{D), as fc —>■ oo. 

We choose an arbitrary w G C^{D), set </> := V • {kVw) and consider the solutions 
fc € N, of the corresponding homogeneous Dirichlet problem (38). Then we have 
for every k G N the trivial identity 

[ dx= [ da; 

J D J D 

and the convergence of the flows follows from the compensated compactness lemma 
[64, Lemma 1.1]. We may thus pass to the limit in the variational formulation to see 
that M = M is the unique solution to the Neumann problem. 

On the other hand, by [54, Lemma 2.2] together with the Holder continuity up to 
the boundary of both k gN, and p, it follows that for fixed x G D, p^^\-,x, ■) 
p{-,x,-) uniformly on compacts in (0,T] x D for all T > 0. It follows from (37) that 
u^^\x) —>■ Uoo{x) for all a; £ H as fe —oo. Hence, u must coincide with Uoa and the 
assertion is proved. □ 

Remark 2. Note that the regularization technique we employed in the proof of The¬ 
orem 6.4 may be easily modified to prove the Feynman-Kac formula 

u{x) = E3,0(Xt-(£,)), X G D 

for the conductivity equation (1) with Dirichlet boundary condition u\dD = </', where 
(j) G and 

t{D) := inf{t > 0 : Xt £ R‘^\D} 

denotes the first exit time from the domain D. Such a proof requires the fact that 

^■ArW(r>) ^•at(d) in law on (^([O, oo); D) 

as n —^ oo for every x G D. This follows immediately from the Lipschitz property of 
dD, implying that all points of dD are regular in the sense of [35, Chapter 4.2], see 
also the proof of Theorem 7.1 in the next chapter. 

A slight modification of the arguments from above yields the following result which 
is in fact a corollary rather to the proof of Theorem 6.4 than to its actual statement. 

Corollary 1. Let f be a bounded Borel function and let a > 0. Then there is a unique 
weak solution u G C(D) n H^{D) to the boundary value problem 

V • (kVu) — au = 0 in D 

subjeet to the boundary condition (3). This solution admits the Feynman-Kac repre¬ 
sentation 

nCO 

(39) u{x) =E^ e-“*/(Xt) dL* for all xGD. 

Jo 

Proof. Repeat the proof of Theorem 6.4, however, substituting {Tt,t > 0} with the 
Feynman-Kac semigroup {Tt,t > 0}, Ttv{x) := Exe~°‘*'v{Xt). Note that in contrast to 
the Neumann problem without the zero-order term, the gauge funetion E^, /g e~* dL^ 
is finite P^j-a.s. for every x G D. □ 


From Feynman-Kac formulae to stochastic homogenization in EIT 


21 


6.2. Complete electrode model. The main result for the complete electrode model 
(1), (5) is the following theorem. 


Theorem 6.5. For given functions f,g defined by (6) and a voltage pattern U £ 
satisfying (1), there is a unigue weak solution u £ C{D)r\H^{D) to the boundary value 
problem (1), (5). This solution admits the Feynman-Kac representation 


(40) 
with 

(41) 


poo 

u{x) eg{t)f{Xt) ALt for all x e D, 

Jo 


eg(t) := exp - y g{Xs)dLs'^, t>0. 


Before we are ready to give a proof of Theorem 6.5, let us introduce the Feynman- 
Kac semigroup of the complete electrode model, i.e., the one-parameter family of op¬ 
erators {Tf,t > 0} defined by 

(42) Tfv{x) := Mxeg{t)v{Xt), x G D and t > 0. 

Let us define the perturbed Dirichlet form ,T>{£^))) by a perturbation of (£,£>{£)) 
with the measure g ■ a, i.e.. 


(43) £^{v,w) :=£{v,w)-\-{gv,w)gjg, v,wGV{£S), 

where 'D{£^) = F[^{D) by the standard trace theorem. 


Proposition 4. The Feynman-Kac semigroup {Tf, t > 0} is a strong Feller semigroup 
on L'^{D). 


Proof. First, it follows from a straightforward computation that {T^,t > 0} is asso¬ 
ciated with the Dirichlet form (£®,T>(£’®)) and is therefore a strongly continuous sub- 
Markovian contraction semigroup on Lf'{D). Tf is a bounded operator from L^{D) to 
L°°{D) for every t > 0, which can be shown using Fatou’s lemma. By the Dunford- 
Pettis theorem, T® can thus be represented as an integral operator for every < > 0, 

(44) T^(j){x)= [ p^{t,x,y)(l){y)dy, for every (j} G L^{D), 

Jd 


where for all t > 0 we have p^{t,x,y) £ L°° (DxD) and p^{t,x,y) > 0 for a.e. x,y G D. 
In order to prove the strong Feller property, we show that Tf,t>0 maps bounded 
Borel functions to C{D). As in the papers [32, 49], we use an iterative method to 
construct the transition kernel density p®: Let p^it, x, y) := pit, x, y) and set 


Pk(J^^^y)'=- / Pis,x,z)g{z)pl_j^(t-s,z,y)da{z)ds, 

Jo JdD 


fc £ N. 


Note that the terms p| are symmetric in the x and y variables by the symmetry of p. 
By induction using Lemma 6.3 it is not difficult to verify that for all £ N 


[ f g{x)pl{s,x,y)da{x)ds < ( sup [ g{Xs)dLs\) 
Jo JdD ^ xGD ^ Jo ’ ^ 


fc+1 


and that there is a positive constant ci such that 

/ r 'i \ fe+1 

(45) |p^(t,x,y)| < sup jEa, / 5(Xs)dLs|j 

^ x^D ^0 


for all fc £ N. 
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Let us show the continuity of /c G N U {0}. For fc = 0, this is a consequence of 
Theorem 1. Now assume that is continuous on {to,T] x D x D toi to > 0, then 
we have for t G {to,T] 



s, z, y) Aa[z) ds 


X, z)g{z)pl_^(t - s, z, y) dai^z) ds. 

'D 

Note that the first integral on the right-hand side tends to zero uniformly as to —0 


which is a consequence of (45), while the second integral is continuous by assumption. 
Hence, there exists a T > 0 such that the series cc, y) := v) converges 

absolutely and uniformly on any compact subset of (0, T] x Dx D and is thus continuous 
on (0,T] X D X D. By the Markov property we have for all t G (0,T] and every x G D 
the following expression for (j){x) 



Therefore, the assertion for arbitrary T > 0 follows from the Chapman-Kolmo-gorov 


□ 


equation. 


Corollary 2. Let u be defined by the Feynman-Kac formula (fO), then u G C{D). 
Proof. Let us define a P^,-martingale by 



where the right-hand side is obtained using the Markov property of X together with 
the fact that eg is a multiplicative functional of X. Obviously, 



Jo 

is a Pa,-martingale as well, and hence we have for all 0 < s < t 


t — S 



eg{s)u{Xs) = eg{s)Ex^eg{t - s)u{Xt-s) + eg(s)Ex, 


eg{r)f{Xr)dLr. 


Setting s = 0 thus yields 



(46) 


for all t > 0. 


Tfu is continuous on D by Proposition 4. To prove that u is continuous on D, it is 
sufficient to show that the second term on the right-hand side of (46) tends to zero 
uniformly in a;, as t > 0. This is, however, clear since we may estimate 



where the right-hand side tends to zero as t —>■ 0 by Lemma 6.3. 


□ 


The following lemma yields a semimartingale decomposition for the composite pro¬ 


cess u{Xt) which compensates for the lack of Ito’s formula in the proof of Theorem 


6.5. 
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Lemma 6.6. Let u € H^{D) denote the weak solution of the boundary value problem 
(1), (5). Then for all t > 0 

f{Xs)dLs+ f g{X,)u{Xs)dLs, 
do 

Px-a.s. for q.e. x € D. 


(47) u{Xt) = u{x) + f Xu{Xs) dM“ - f 
Jo Jo 


Proof. Applying the Fukushima decomposition (26) to the perturbed Dirichlet form 
(f®, H^{D)), we obtain the unique decomposition 

v{Xf) - v{Xl>) = Mf’" + Nf’^, for all t > 0, P^-a.s. for q.e. x GD 


into a martingale additive functional of finite energy and a continuous additive func¬ 
tional of zero energy of the non-conservative Hunt process X^ associated with {£^,H^{D)). 
We study the relation between the continuous additive functionals N'" and N^’’". Let 
us assume first that v is in the range of the 1-resolvent associated with the perturbed 
Dirichlet form, i.e., 

v{x) = Gl(l){x) = Ex 

for some (p G Lf{D). Then we have by the resolvent property, —CG\(j) = (f> — v so that 
for all w G H^{D) 

£^{G\(j),w)= f {cj) — v)wdx. 

Jd 

Using the Revuz correspondence, we see that iV®’’' admits a semimartingale decompo¬ 
sition, namely 

Nr= A</>(Xf)-u(Af))ds. 

^0 

Moreover, an easy computation yields that 

(48) Gi(l>-G^^(j) = Ex{J^ e“‘5(At)G?(^(At)dLt} for all a > 0. 

By Corollary 1, the right-hand side in (48) is the unique weak solution of the elliptic 
boundary value problem 




X G D 


I V • {kVw) — w = 0 in D 

= gG\(j) on dD. 

That is, it coincides with the 1-potential Ui{Gl(j){g ■ a)) G H^{D). Invoking the Revuz 
correspondence once more, we see that the zero energy continuous additive functional 
in the Fukushima decomposition of the 1-potential corresponding to the signed Radon 
measure vg ■ a is given by 

rt pt 

t/i(ug-a)(A,)ds- / u(X«) 5(A:,) dL„ for alH > 0. 

Jo Jo 

Therefore, we obtain 


Nf = TVf - W 


j-Ui{vg-a) 


f 


{(j){Xs) - Gi(j){Xs))ds - / Ui{vg ■ <T){Xs)ds + / v{Xs)g{Xs) dLs 
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Moreover, notice that is related to X by a random time change, namely 


X? = 



s < C® 
s>C®, 


where the lifetime C® is given by 

C® := inf : J g(Xs) dLs > 

and Z is an exponentially distributed random variable with parameter 1. Hence, we 
obtain 


(49) N- = Xr + f v{XI)g{XI) dL, for all t < C®. 

do 

This equality may be generalized to the case of an arbitrary v € H^{D) not necessarily 
in the range of the resolvent using an approximation argument. Namely, we consider the 
sequence (ufc)fegN with Vk '■= = G\<l)k^ (f'k ■= fc(u — Then Vk £ H^{D) 

for all fc G N and the sequence {vk)kGN satisfies both. 


lim £^{vk — v,Vk — v) = 0 and lim £(vk — v,Vk — v) = 0 

k—>oo k—^oo 

so that by [25, Corollary 5.2.1], there exists a subsequence, for convenience still denoted 
(i’fe)fceN, such that Vk{Xf) —>• v{Xf), —>• iVf’’' and —>• N]’ uniformly on each 

Hnite time interval, Pa;-a.s. for q.e. x € D. In particular, it follows that (49) holds for 
arbitrary v € H^{D). 

As u solves the boundary value problem (1), (5) we have that 

£’®(u, v) = (/, v)gjj for all v £ H^{D) fl C{D) 

Since the perturbed Dirichlet form {£^,H^{D)) is regular, we obtain from the Revuz 
correspondence the representation 

A^f’“ = — f f{X^)dLs, Pa;-a.s. for q.e. x G D. 

Jo 

Finally, using the representation 
M“ = 

for the martingale additive functional, the claim follows from the Fukushima decom¬ 
position (26), (49) and the Markov property of X. □ 


Vu(Xs)dM“, Pa;-a.s. for q.e. x G D 


Proof of Theorem 6.5. There exists a weak solution u G H^{D) of the boundary value 
problem (1), (5) so that with regard to Corollary 2, it remains to show that this weak 
solution u admits the Feynman-Kac representation (40). Note first that the gauge 
function 

POO 

(50) E, / egit)dLt 

Jo 


is finite P^j-a.s. for every x G D, hence the expression on the right-hand side of (40) is 
well-defined. Lemma 6.6 yields the semimartingale decomposition 

uiXt) = u{x) + [ VuiXs) dM“ - 
Jo 


fiXs) dLs 


g{Xs)u{Xs)dLs 
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Pa^-a.s. for q.e. x £ D. Note that the second term on the right-hand side is a local P^;- 
martingale and that Cg is continuous, adapted to > 0} and of bounded variation. 

Multiplication by such functions leaves the class of local martingales invariant. Using 
integration by parts we thus obtain for q.e. x £ D the identity 

u{Xt)eg{1:.) = u[x) + f eg(s)Vu(Xs) dM“ - f eg{s)f{Xs)dLs, 

Jo Jo 

Px-a.s., where the second summand on the right-hand side is a local P2;-martingale. 
That is, there exists an increasing sequence of stopping times which tend to 

infinity such that for every k £ N 

ptATk 

MtAr,-= eg{s)Xu{Xs)dM: 


is a P^j-martingale. By definition of the term eg, it is, however, clear that 
Ej; sup|Alt/\T-^| < oo for alH > 0 and every x £ D 

fcGN 

which is sufficient for {Xit ,t> 0} to be a Pa;-martingale by the dominated convergence 
theorem. Hence, 

u{x) = Ea; / eg{s)f{Xs)dLs +E,^u{Xt)eg{t) for q.e. x £ D. 


Letting t ^ oo finally yields 

poo 

u(a;) = Ex / eg{t)f{Xt)dLt for q.e. x G £>, 

Jo 

where we have used the fact that u is essentially bounded by standard elliptic regularity 
theory. Finally, by the fact that we have actually shown in Corollary 2 that the right- 
hand side in the last equality is continuous up to the boundary, the assertion holds for 
every x £ D. □ 

Remark 3. Note that the technique we used to prove Theorem 6.5 fails for the Neu¬ 
mann problem corresponding to the continuum model. This comes from the fact that 
in this case the gauge function (50) becomes infinite. For the same reason Theorem 1.2 
from [14], specialized to a zero lower-order term, does not yield the desired Feynman- 
Kac formula for the continuum model either. 

6.3. Mixed boundary value problems. Now we can directly deduce the desired 
Feynman-Kac formula for the mixed boundary value problem corresponding to the 
stochastic anomaly detection problem introduced in Section 3.2 of the previous chapter. 
Recall that in this setting dD consists of two disjoint parts diD and 82D and that 
measurements can be taken only on the accessible boundary diD while the electric 
potential vanishes on the inaccessible boundaryd2D. The deterministic FIT forward 
problem for the complete electrode model is then given by the conductivity equation 
(1) subject to the mixed boundary conditions 


(51) 


K.V ■Vu\dD + gu\dD = f ondiD 
u\dD = 0 on 820 . 


The following result is a corollary to the line of arguments that led to the proof of 
Theorem 6.5 rather than to its actual statement. 
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Corollary 3. For given functions f,g defined by (6) and a voltage pattern U G 
satisfying (7), there is a unigue weak solution u € C{D)C\Fl^{D) to the boundary value 
problem (1), (51). This solution admits the Feynman-Kac representation 

(52) u{x) = f eg{t)f{Xt) dLt for all x G D, 

Jo 

where r := inf{t > 0 : Xt G 02 -D} denotes the first hitting time of 820. 

Proof. Repeat the computations from Subsection 6.2 with the Feynman-Kac semigroup 
{Tf,t > 0}, where Tfv{x) := IEa;{[i < deg(t)u(Xt)} instead of {Tf,t > 0}. □ 

7. Stochastic homogenization 

In this section, we show that the stochastic FIT forward problem may be homoge¬ 
nized both theoretically and numerically by homogenization of the underlying diffusion 
process on the whole space Moreover, we provide a continuum version of a quan¬ 
titative estimate which has been obtained recently for the discrete random walk in 
random environment in [27, 19]. 

7.1. Preliminaries. For convenience of the reader, let us recall some standard con¬ 
cepts from homogenization theory. Let ■= G i = l,...,d, 

denote a vector field. We say that <)> is a gradient field if for every if G 

/ (fidjtp — dx = 0, i,j = l,...,d. 

Jr’^ 

Moreover, we say that (j) is divergence-free if for every if G 

(fidifi da: = 0. 

Now let us consider a conductivity random field {k{x,u!), {x,u!) G x F} and let 
{©X ,x G R"^} denote the underlying dynamical system which is assumed to satisfy 
the assumptions (i)-(iv) from Subsection 3.2. A vector field (f> G L^(r;R‘^) is called a 
gradient field, respectively divergence-free, if its realizations (j){-,uj) : 1 —>■ x 1 —>■ 

are gradient fields, respectively divergence-free, for P-a.e. w G F. We define 
the function spaces 

-^pot(r) := {</> G L^(F;R‘^) : is a gradient field P-a.s.} 

LgQi(r) := {0 G L^(r;R‘^) : </>(•, w) is divergence-free P-a.s.}. 

li <f> G Lpot(F), we can find a function 77 : R"^ x F —>■ R such that 'g{-,Lo) G and 

(53) Vr]{-,uj) = a.e. in R'^ for V-a.e. w G F. 

In particular, (53) defines a stationary random field with respect to the measure V. 

We call 77 the potential corresponding to 4>. 

Remark 4. Note that 4> G Lp^^(r) does not imply that {(/>(x,uj), (x,uj) G R'^ x F} is 
a stationary random field with respect to P. In fact, it can be shown that this is not 
true for d = 1. 
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We define another function space 

Vl,, := {<A e : M0 = 0}, 

so that one obtains an orthogonal Weyl decomposition of namely 

cf., e.g., [64]. Let ^ denote a direction vector, i.e., j^j = 1. The so-called 

auxiliary problem for the direction ^ reads as follows: Find € VpQt(r) such that 
G -^Li(r) or equivalently, 

(54) .</,} = 0 for all cj) € Vpot(r). 

For a proof of existence and uniqueness of the solution to the auxiliary problem we 
refer the reader to the seminal paper [48] by Papanicolaou and Varadhan. 

We can now bring the underlying diffusion processes evolving in the random medium 
into play by recalling a stochastic homogenization result which was obtained by Lejay 
[42]. Let {Ke{x,uj), {x,uj) G x F} denote the scaled conductivity random field, see 
Subsection 3.2, and let denote the diffusion process on which is associated 
with the regular symmetric Dirichlet form 


£‘^’^v,w) 



Ke(-, w)Vz; • Vw dx. 


v,we V{£^'^) := 


on It has been shown in [42] that, under assumption (Al) from Subsection 

3.2, for V-a.e. w G F 

(55) —>■ X* in law on C'([0,oo);R‘^) as e —>■ 0, 

where X* denotes the diffusion process on which is associated with the homogenized 
Dirichlet form 


£*{v,w):= [ K*VvVwdx, v,w&V{£*) ■= 

Js.<‘ 

on L‘^(W^) and the constant, symmetric and positive definite matrix k* satisfies the 
equation 

(56) C-K*C = M{(^-hx^)-'«(C + X^)}, 

where G Vpot(r) denotes the solution to the auxiliary problem (54) for the direction 
^ G 

7.2. Homogenization of the EIT forward problem. The following theorem is 
our main result. Its assertion is in fact a rather direct consequence of an invariance 
principle for reflecting diffusion processes obtained recently by Chen, Croydon and 
Kumagai [13] and the Feynman-Kac formula (52) from Corollary 3. 
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Theorem 7.1. Let {Ke(a;,a;), {x,ui) S R'* x T} be a stationary random field satisfying 
assumption (Al) from Subsection 3.2 and assume that the trajectories satisfy k{-,uj) € 
or piecewise constant for V-a.e. u G F; let S = 0. Then, for a 

given voltage pattern U G R^, we have for the potentials in the stochastic boundary 
value problem (11), (12) 

(57) Ue{x,iij) —>• u*{x), X £ D for V-a.e. a; G F, as e — >■ 0, 
and the corresponding electrode currents satisfy 

(58) lim J/(e, w) = r—r f k*^ ■ \7u*{x)\g.D dcrlx) for V-a.e. w G F, 

^^0 \Ei\ 

I = 1, ...,N, where the function u* G Hq{D U diD) fl C{D) is the unique solution to 
the deterministic forward problem 

(59) V • (k*Vm*) =0 inD 
subject to the boundary conditions 

■yu*\g^D + gu*\diD = f ondiD 

u*\d 2 D = ^ on d 2 D 

with a constant, symmetric and positive definite matrix k* given by (56). 


Proof. Let us first show that for V-a.e. w G F, Ue{x,uj) —>■ u*{x), a: G D, as £ —>■ 0. Let 
(e^fc)fcgN be an arbitrary monotone decreasing null sequence and let X‘^’^ denote the 
reflecting diffusion process on the half-space corresponding to the regular symmetric 
Dirichlet form iJ^(Rl U R'^"^)) on L^(Rl UR"^”^). By assumption (Al) from 

Subsection 3.2, we deduce from [13, Section 4] that for V-a.e. w G F 

^u,ek Qg C([0, oo); R1 U R^^”^), as fc — >■ oo, 

where X* is the reflecting diffusion process on the half-space associated with the ho¬ 
mogenized regular symmetric Dirichlet form {£*, H^{Wf UR‘^“^)) on L^(Rl UR‘^“^). 
The constant, symmetric and positive definite matrix k* is given by (56). 

Let us first show that for every x € D and V-a.e. realization w G F of the random 
medium 

X^^t.ek ill on C'([0,oo);R1 UR^^”^), as fc —>• oo. 

Consider the functional F : C'([0,oo);Rl UR‘^“^) —>■ [0,oo], 


(j) 


oo, 

inf{t > 0 : |()'(t)| 


R} 


if for alH > 0 : \(j){t)\ < R, 
else, 


defined in such a way that F{X‘^’^>=) = Let {(j)k)ke« denote a sequence of 

continuous functions that converges uniformly towards </> on compacts in [0, oo). If 
liminffc_>.oo T'((/)fc) is finite, then we may extract a subsequence, still denoted {(j)k)keN 
for convenience, such that F{(j)k) —>■ liminffc^oo F{()k)- We have 

|())(lim inf F((;ifc)) - (^fc(F(())/c))| <|())(lim inf F{(l)k)) - 4>iF{(f>k))\ 

k—¥<x> k—¥oo 

+ \cj{F{cjk))-MF{<(k))\ 


and by our assumption, the right-hand side vanishes as fc —^ oo. From the closedness 
of diD, we conclude hence that ^(liminffe_>oo A(^fc)) G diD. In particular we have 
that F{(j)) < liminffc_>oo F^ejk) if sup^^jp-j.]|^(t) — (t)k{t)\ —0, as t —>■ oo, for all T > 0, 
i.e., F is lower semi-continuous. Now assume that is a discontinuity point of F, i.e.. 
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there is <5 > 0 and ko such that F{(j)) + i5 < F{cj)k) < oo for all k > ko. Then it follows 
that G D for all t G liminffe_j,oo However, the boundary diD is 

regular in the sense of [35, Chapter 4.2], that is, a diffusion process originating from 
X G d\D will immediately exit from D P“’^''-a.s., respectively P*-a.s. In other words, 
the set of discontinuities of f is a null set with respect to these measures and hence 
the claim follows from the continuous mapping theorem, cf. [8]. Now, as in the proof 
of [53, Theorem 5.1], it follows with the Fukushima decompositions in Section 7 of the 
second chapter that 




Ar^’ 


T^,^k 




in law on (^([O, oo); Ml U ^ x R_|_), 


as fc —oo. With regard to the Feynman-Kac formula (52), the assertion (57) can be 
proved as follows: 

Suppose first that / is continuous. Then we have to show that 


lim sup 

j-s-oo 





-E, / egit)f{Xt)dLt 

Jo 


= 0 . 


We estimate this with a difference of truncated Riemann sums 


I ^ (e, / e-’^A4)/(^r/0[4+i <r,]dLr^ 

fc=0 

f eg{tk)f{Xti^)[tk+i<T]dLt), 

Jtk ' 

where tk = kh and h is a step size. The difference of truncated Riemann sums for 
fixed X and h goes to zero as j —>■ oo if we assume that ) converge 

weakly to {X,L,A), where = log egand A = log eg. 

The error terms for and X are nearly analogous, so it is enough to consider 

just in detail. The increments (excluding the edge case where tk < tj < tk+i 

which we will omit but which can be treated in the same way) are of following form 



e-’^^(t)f(Xr^) - e-’^A4)/(Xr/^)) dL'. 


i^,ej 


fik + l 

/ {e-g’^^{t)-e-’^^{tk))f{XrndL' 

Jtk 


+ eg’"A4) dL^’^^ 

tk 

The latter term can be handled with the continuity of the paths of X together with 
the uniform continuity of / on the compact set D. This is since 


ptk+i 

eg’^AA) / (/(Xr"0-/(^rO)dLr’"^ 

Jtk 

< ( 2(1 - M0h{x-'^n))\\f\\^ + 02s{f)MOh{x‘^’^n)) 

where 0 s(x) is the maximum variation of the function x on the interval [0,X] 

0s{x) := sup { \x{t) — x(s)l; 0 < t, s < X, jt — s] < (5} 

and tps (t) is the continuous approximation of the indicator function [t < 5] with support 
in [—26,26]. Therefore, after taking the limit j —>■ 00 , the latter terms give that the 
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corresponding total approximation error can be bounded by 

4\\f\\^E4[9hiX) > S]L^^N} + 2e2s{m,{LrAN} 

which goes to zero as <5 —>■ 0. 

The first term can be estimated by 

rtk+i 

•ttfc 

ftk+i 

< 2ii/iiooii5iioo / - L^n = \\f\u\ 9 \uLZ[ - zn 

Jtk 

since g > 0. Therefore, after taking j —>■ oo the first term gives the total approximation 
error that is bounded by 

2II/II00II5II00 ^Ea;{[0;^(L) > S]LtAn} + 2 ( 5 E 3 ,{Lt-/\ 7 v}^ 

which goes to zero as <5 —>■ 0 as well. 

We are still left with the truncation. The truncation can be removed since Tj and r 
are a.s. finite and moreover, converges weakly to and thus, we get a uniform 

estimate 


limsup Ex 

j->oo Jo 

< ‘2\\f\\ooEx{Lr — LxAn} 


eZ^it)fiXZ^)dLZ^ - eg{t)fiXt)dLt 


which gives the claimed convergence for continuous / if we can prove weak convergence 
of the sequence 

Therefore, we will next verify the weak convergence of jointly with 
We can assume that X and —>• L almost surely in (7(0, T). Suppose that 

g = 'ZifLi where gi > 0 and continuous, Eif\Eu =% and cr{dEi) = 0. Therefore, 
it is enough to show the convergence for g = [Ei\ for Ei open in dD. Since Ei is open, 
it can be approximated from below by an increasing sequence of continuous functions 
g'j^ that converge pointwise to g. 

It follows from the almost sure convergence of and that 

rt pt pt 


liminf / > limsup / g'^.{Xt)dLt= / g{Xt)dLt. 

j^oc Jq fc-j-oo Jo Jo 

Since and X will not hit dEi in [0, T] almost surely, we can get the other direction 

by considering 1 — 5 . Therefore, 

pt pt 


limsup/ g{X^'''^)dLt’‘'^< g{Xt)dLt. 
j-tco Jo Jo 


Using countability of the rational numbers, we obtain 


Vt e [0, T] n Q: lim {t) = A{t) 

j->oo 

almost surely and by monotonicity this implies the almost sure convergence of to 
A in (7(0, T), which implies the weak convergence of the original versions. The same 
technique allows us to extend the assertion to the case of discontinuous /. 

For the proof of (58) note that the boundary condition (12) allows us to write 

(61) Ji{ek,u;) = f if - gUeki-,Z\diD)da{x) for P-a.e. w £ F 




From Feynman-Kac formulae to stochastic homogenization in EIT 


31 


and that (12) may be written in the form +gl)u^ = f. We deduce from the well- 
posedness of the forward problem that L‘^{dD) —>• L'^{dD), is bounded. 

Since, due to assumption (Al), the corresponding constant does not depend on e, 
the sequence {ue^)keti is bounded in LF'{dD) for P-a.e. w G F implying its uniform 
integrability, see, e.g., [41]. As we already know the pointwise convergence Uef.{x,uj) —>■ 
u*{x), X € D, as k ^ oo, an application of Egorov’s theorem yields convergence in 
L^{dD) so that the assertion follows by the triangle inequality and taking the limit 
k ^ oo inside the integration in (61). □ 


Remark 5. We would like to point out that the effective conductivity k* is determined 
by the invariance principle on the whole space 


Remark 6. We would like to mention that for more general functions / we can prove 
that the well-posedness of the forward problem (11), (12) together with the weak 
convergence of A“’l) and some stopping arguments give a slightly weaker 

form of (57), namely that the convergence holds in P-probability. 


7.3. Continuum approximation of the effective conductivity. In this subsec¬ 
tion, we provide the theoretical foundation for the convergence analysis of numerical 
homogenization methods based on simulation of the underlying diffusion process. More 
precisely, a rigorous convergence analysis of such a method requires a quantitative esti¬ 
mate that is stronger than the qualitative result (55), which was obtained in [42] using 
merely the central limit theorem for martingales. We provide such a quantitative re¬ 
sult in the following theorem by generalizing a classical argument due to Kipnis and 
Varadhan [40]. The proof relies on new spectral bounds bounds which were obtained 
recently by Gloria and Otto [29]. We refer the reader to the recent papers [19, 27] 
for an analogous estimate for the discrete lattice random walk in random environment 
as well as to the paper [44] which was the first one to use the Kipnis and Varadhan 
argument in order to obtain quantitative results. 


Theorem 7.2. Let {ki:{x,u!), {x,u!) G ^ T} be a stationary random field satisfying 
assumptions (Al) and (A2) from Subsection 3.2. Then for every direction vector G 
R"^ there exist positive constants ci, C 2 such that 


(62) 


2 t 




< Cl 


|logt|°^t ^ 


d = 2 
d= 3. 


Proof. For hxed w G F let us consider the diffusion process X“ on R'^ which is asso¬ 
ciated with the symmetric regular Dirichlet form {£‘^d^ on L^(R'^) under the 

measure Pg . Following [40], we search for a decomposition of the form 

(63) Xf = Mr + Rf, 

where Mf is a Pj)"-martingale and for every direction ^ G R'^ the projected remainder 
Rf ■ f converges to zero in A^(F) as t —^ oo. 

Once we have found a suitable decomposition (63), we first show that 

(64) - {M, ■ on = t-m{R, ■ 0^ 

Then we will use spectral calculus to estimate the right-hand side of (64). 

In order to obtain the decomposition (63), we recall that the auxiliary problem (54) 
is equivalent to the following stochastic elliptic equation in the physical space R"^: Find 
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G Vpot(r) such that for V-a.e. w G F, the corresponding potential 77 ^ : x F —R 

is in C(M‘^) fl as a function of x, satisfies 77 ^( 0 , w) = 0 and 

(65) - V • K(-,a;)(C +V77«(-,a;)) = 0 in R'^. 

Let us define the function 


()) : R'^ X F —>• R'^, (j){x, w) := x + r]{x, w) — 77 ( 0 , w), 

where 77 := ( 77 ®^ and r]^\ i = denotes the potential corresponding to 

the solution to the auxiliary problems for the coordinate directions. As the transition 
density kernel of X‘^ is jointly Holder-continuous and G C'(R‘^) fl 

7 = 1 , d, for V-a.e. u; G F, the Fukushima decomposition of (j)i{X^,uj) holds for every 
starting point x G R'^ rather than quasi-every x G R"^. A straightforward computation 
using the fact that 77 is defined via the auxiliary problem yields that the continuous 
additive functional of zero energy in the Fukushima decomposition vanishes so that 

= 7 = i,...,d. 

We set 

Mr and -- 77 (Ar, 0 .) + 77 ( 0 ,cc) 

and consider the quantity 

(66) l(w • O' = E(M* • O' + nR't ■ o' + 2l(M, • Oi.Rt ■ 0- 

By computing the predictable quadratic variation of the martingale additive functional 
we obtain for alH > 0 and a.e. x G R"^ 

Er (Mr • o' = Ex ( 2(e + V77«(A„ cc)) ■ 77(A„..)(C + V77«(A„..)) ds). 

By the stationarity of {V? 7 ^(x,u;), {x,u}) G R"^ x F} with respect to V we have thus 
I(Mr • O' = 2 ^ • for all t > 0 . 

Moreover, as in [18], it follows that the last term on the right-hand side of ( 66 ) vanishes. 
This can be seen by studying the so-called environment seen by the particle process, 
i.e., the stochastic process defined by 

‘la; t = 0. 


{Y^^, t > 0} is a stationary process with respect to the annealed measure P, i.e., for ev¬ 
ery finite collection of times t^'^\ * = 1, k, the joint distribution of y)(i)_|_^,..., y((fc)_|_^ 
under P does not depend on h > 0. It is well known that the underlying dynamical 
system {©X , X G R'^} defines a d-parameter group {Sx,x G R'^} of unitary operators 
on L^(F) by '■= '0(©xw) and this group is strongly continuous, cf. [42]. Its d 

infinitesimal generators {Di,'D{Di)), ..., {Dd,'D{Dd)) are given by 




lim 
h —^ 0 -|- 


Sheilp - -0 
h 


i = 1, ...,d, 


for all 0 G L^(F) such that the limit exists. These operators are closed and densely 
defined. We denote D := {Di,..., D^)’^ and introduce the infinitesimal generator 
{C,'D{C)) on L^(F) of the environment viewed by the particle process, that is, the 
non-negative definite self-adjoint operator C := —D ■ kD on L^(F). Note that due 
to the fact that the trajectories of the random conductivity field are in general not 
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differentiable, one has to use Dirichlet form theory in order to give a precise meaning 
to computations involving C. We refer the reader to the work [42] where this has been 
carried out in some detail. By the self-adjoinedness of C,^ the law of the environment as 
viewed from the particle process under P is invariant with respect to time reversal and 
is odd by [21, Corollary 2.1], i.e., it changes its sign under time reversal, whereas 
i?“, which is the zero energy part of the Fukushima decomposition (63), is even by [21, 
Theorem 2.1]. Thus, we notice that the identity (64) holds, as claimed. 

We will next show that the estimate (62) follows from the following spectral repre¬ 
sentation 


(67) 


E(i?t • O' = 2 / (1 - e-^‘)A-2 u«). 


where {Ex^ A € K} is the unique spectral family given by the spectral theorem such that 
C = XdEx and the function := D ■ G Indeed, given the formula (67) 

for the projected remainder and due to the assumption (A2), we can now exploit the 
following optimal estimate from [28, 29]: For all 0 < 7 < 1, there exists a positive 
constant c such that 


f d{Exv^,v^) < 

Jo 


More precisely, we split the integral (67) into three parts, the first ranging from 0 
to the second from to 1 and the third from 1 to 00 , respectively when t > 1. 
For the latter we have the trivial estimate 

/ OO pOO 

d{Exv^,v^)< J d(i5Ai’^, 1’^) = M{(d^)^}. 

The first part is bounded by a positive constant as well, namely by 

poo 

/ d{Exv^ = t / / 

Jo Jo J\ 


= t 


poo pocAt 

/ a~‘^ / d{Exv^,v^) 

Jo Jo 

poo 

< ct a~^{a da. 

Jo 

Similarly, the second part can be estimated by 

pi pi poo 

/ X~^ d{Exv^,v^) = 2 / q;“^ dad(i5Ai’^,'i’^) 

7t-i Jt-i J\ 

poo 

< 2c a“^(a A 1)^*/^^^ da, 

Jt-^ 


which diverges logarithmically for d = 2 and is bounded by a positive constant for 
d = 3. Therefore, combining these computations with the identities (64) and (67) the 
estimate (62) follows. 

It remains to prove the spectral representation (67). In order to take advantage of 
the spectral theorem, we would like to express the projected remainder in the form 
for some bounded continuous functions ipi and ^/' 2 . However, for 
this we would need the components of rj to be stationary, which is not the case. Fur¬ 
thermore, we would want to use ipiY) — which is unbounded at zero. 
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Inspired by the computations in [48], we can find a remedy for both of these ob¬ 
structions; namely we consider the function 

:=-?75(X“,w)775(0,0;), 

where r]s is defined in analogy to rj with the difference that it corresponds to a different 
auxiliary problem, modihed by a zero order term. This modified auxiliary problem 
reads as follows: Find 77 ^(-, 0 ;) G C'(IR‘^) fl ^ > 0 such that for V-a.e. a; G F, 

the random field {r]g{x,uj), {x,uj) G x F} is stationary with respect to V with 
Mr]g{x, •) = 0 for every x G R'^ and satisfies 

for V-a.e. w G F. We refer to [48] for the proof of existence and uniqueness of 77 ^, 
(5 > 0. Note that 77 ^( 0 , w) 7 ^ 0 in general so that we have for the projected modihed 
remainder the expression 

(68) E{K/ ■ 0^ = VsiXt)? - 2l77«(X„.)77|(0,-) +l(77f(0,-))". 

The equivalent formulation on T^(r) of the modihed auxiliary problem for the di¬ 
rection ^ G R'^ reads as follows: Find the unique solution 77 ^ G T^(F) of the elliptic 
equation 

Sr/l — D ■ -|- Drjg) = 0 in F. 

In particular, the function was chosen such that = (i5 -I- C)ril = Crj^. 

We will now only need to show that the modihed remainder for hxed 6 and hxed t 
can be written in the form M.'ipi{C)v^ip 2 {C)v^, where 7 /’i(x) = 2{x-\-d)~^ and '<p 2 {x) = 
2e“‘“(x -I- 5)“^. In fact, the introduction of the zero order perturbation removes the 
singularity coming from the term x~^. 

The hrst term on the right-hand side of ( 68 ) may equivalently be written as 

1(77«(X„.))' =1(771(0, y,))2 =M(77^(0,•))^ 

where we have used the stationarity of the environment as viewed from the particle 
with respect to P and the stationarity of the random held {r]^{x,uj), (x,uj) G R'^ x F} 
with respect to V, respectively. Therefore, treating the second term on the right-hand 
side of ( 68 ) we obtain 

(69) l(i?,-' • = 2M(77|(0,-))" - 2M77«(0,-)T;'775«(0,-), 

where > 0 } denotes the strongly continuous semigroup on L^(R^) associated 

with i7^(R'^)) which satishes 

E^rjl(X,-,co)=Trr,l(0,L,). 

Going back to (69), respectively the corresponding identity on T^(F), we have thus 
for the hrst term on the right-hand side 

pOO 

2 M(r 7«)2 = 2M(5 -h £)~^7;^(5 -t- = 2 ((5 -h A)-^ d{Exv ^,««), 

Jo 

whereas the second term may be written in the form 

poo 

2M(5 + C)~^v^e~^^{6 + C)~^v^ d{Exv^,v^). 

Jo 
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Altogether we have obtained 

POO 

E{r/ -0^ = 2 (1 - e-^*){S + A)-2 d{E^vi, v^) 

JQ 

and the right-hand side of this equality converges as d —0 so that the spectral repre¬ 
sentation (67) holds. □ 

We would like to remark that it appears that the spectral properties of operators 
with stationary coefficients and the stationarity of the corresponding fields (as in the 
proof of Theorem 7.2) might be very closely related in general. However, analysing 
this potential connection that is touched by the work [48] is left for the future studies. 

8. Conclusion 

We have derived Feynman-Kac formulae for the forward problem of electrical im¬ 
pedance tomography and studied the interconnection between these formulae and sto¬ 
chastic homogenization. Using the properties of the underlying diffusion processes and 
some new spectral estimates from [27, 28] we have then obtained a bound on the speed 
of convergence of the projected mean-square displacement of the processes. These re¬ 
sults provide the theoretical foundation for the development of new scalable continuum 
Monte Carlo homogenization schemes. 

Both, the homogenization of the forward model for the complete electrode model 
and the stochastic numerical approximation of the effective conductivity have direct 
applications in EIT anomaly detection problems for random heterogeneous background 
media, cf. [55]. 
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